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1. Introduction 



Shape invariance has proven to be a useful paradigm for understanding exact solv- 
ability in quantum mechanics || ||. It is of interest to understand how broadly shape 
invariance can be applied, and in this paper, I will identify new arenas in which it can be 
used effectively. 

In the past few years, there has been some consideration of extensions of quantum 
mechanics designed to comport well with string theory |J. In particular, string theory 
predicts that there can be a modification in the familiar Heisenberg uncertainty relationship 
||], so that it takes the form 

Ax>-^- + bAp . (1.1) 
Ap 

This relationship is sometimes termed a "minimal length uncertainty relationship," since 
it implies a lower bound for Ax, namely Ax > y a/b. Of course, in order for such a 
relation to appear in the non-relativistic limit, the Schrodinger equation must be modified. 
It is precisely such modifications that are considered in 0, in which some exactly solvable 
models are found that satisfy the minimal length uncertain relationship ( |1 . 1| ) . 

It is the point of this letter to show that the exact solvability of these models is no 
accident. These models turn out to be representatives of a new class of shape invariant 
models. In this letter, I will identify these new shape invariant models, and use the shape 
invariance to obtain solutions. One will see that these models include the models considered 
by @ , as well as some spectral problems that generalize the Schrodinger equation in ways 
beyond the modification ( |1 . 1| ) ■ An explicit consequence of this result is to re-affirm the 
close connection between shape invariance and exact solvability for spectral problems in 
general, and not simply in the context of the Schrodinger equation. We also make technical 
progress, demonstrating that shape invariance can be applied productively in a setting in 
which second-order differential operators can be written as the product of pairs of first- 
order differential operators that are not adjoints of each other. 

2. One-dimensional generalized harmonic oscillator 

We begin by displaying the simplest model of 0. To begin, H considers a modified 
commutation relation 

[x,p] =ih(l+f3p 2 ) . (2.1) 
1 



Such a commutation relation leads to the uncertainty relationship 

h/ 1 



A ^2^ + ?AP J ' (2 ' 2) 
This commutation relation can be realized in the momentum representation by operators 

d 



x = ih 



p = p 



+ 1P , . 

(2.3) 



We note that the parameter 7 is a trivial addition here; its presence does not serve to 
modify the conventional commutation relation, and it is easily set to zero by a canonical 
transformation. It is included here not simply for making contact with [Q, but more 
importantly to lay the groundwork for one of the considerations needed for the multi- 
dimensional oscillator considered in the next section. 

Quantizing the classical one-dimensional harmonic oscillator 

H =^P 2 + \^ 2 x 2 (2.4) 

according to (|27l]) produces the momentum space Schrodinger-like time-independent equa- 
tion 

2E 



»hu{((l+Pp 2 )— ) +2 7P ((l+[3p 2 )—)+ 7 ([3+ 1 )p 2 + 7 } + —p 2 ]y(p) = — 

(2-5) 



Remember, too, that the physics cannot depend on 7. 

We now wish to demonstrate that this equation can be solved exactly via shape 
invariance techniques. This provides a simpler understanding of the exact solvability (and 
a simpler calculation of the energy levels and wavefunctions) demonstrated by analytical 
means in 0. 

A quick review of shape invariance is in order here; more detailed treatments are 
found in the literature, including a review || and a recent article revealing the underlying 
algebraic structure of shape invariance as a BPS phenomenon associated with centrally 
extended supersymmetry |J. 

Suppose one has a Hamiltonian H = AA for which the spectrum is bounded from 
below by zero. In the familiar applications, A = A*, which enforces the positive semi- 
definiteness of the spectrum, but as long as the spectrum is positive semi-definite, A need 
not be the adjoint of A. The eigenfunctions for such Hamiltonians with eigenvalue zero are 
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those annihilated by A. The partner Hamiltonian H = A A will have the same spectrum 
of non-zero eigenvalues as H, since if H(p = Kef), then H(Acf)) = A(A<p). Typically, Atfi = 
has solutions and A<p = does not, and the spectral degeneracy does not extend to the 
zero eigenvalues. 

Shape invariance exists when H and H have the same mathematical form, so that 
H(g) = H(g') + c, where g denotes the coupling constants of the theory, g' are the trans- 
formed coupling constants which are functions of the g, and c is some c-number, which also 
generally depends on the parameters of the theory. This means that the energies of the 
two Hamiltonians not only exhibit the degeneracy referred to in the preceding paragraph 
(Ek+i — Ek), but also satisfy the relation Ek(g) = Ek(g') + c. These results enable us to 
find the energy eigenvalues algebraically. 

The wavefunctions are also found algebraically, through standard raising and lowering 
operator techniques. The ground state of H is the solution of A<p = 0; the states <pk of 
H are related to the states <pk of H by 4>k(g) = <t>k{g') as well as by <fik+i = A<fik, which 
allows the eigenstates to be built up algebraically from the ground state of H. 

We now demonstrate the solution of ( |2.5| ) via shape invariance. To begin, define the 
operators 

d 

A = (a + bp 2 )— + cp 

* (2.6) 

A = — (a + bp )— — h cp . 
dp 

The real numbers a, 6, and c are all positive. Note that with this convention, there is a 
normalizable state annihilated by A but not one annihilated by A. 
Consider the equation TC(p = A0 where Ji = A A. Then 



H(a, 6, c) 
Similarly, 7i = AA is 

7i(a, 6, c) 



(a + bpij- 



cT 2 



+ c(c-b)p 2 -ac . (2.7) 



{a + bp 2 )^- 
dp 



cT 2 



+ c{c + b)p 2 + ca . (2.i 



Consequently, 7i and 7i form a shape invariant pair, since 

H(a,b,c) = H(a,b,c + b) + 2ca + ba (2.9) 
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Applying the principles of shape invariance, one readily obtains the eigenvalues for 
H(a,b,c)4> k (p) = Afc<Mp) and Tt{a,b,c)<p k (p) = A&0&(p). One finds A = 0, A = 
a(b + 2c) = Ai, Ai = a(b + 2(c + b)) + a(b + 2c) = A 2 , and so forth, resulting in the 
general expression 



A fc = a(k 2 b + 2kc) 



k = 0,1,2,3,... 



;2.io) 



We can easily make contact with previous results. Note that, as pointed out above, the 
parameter 7 has no physical content, and thus may readily be set to zero. Once this is done, 
the modified harmonic oscillator equation (|2.5| ) of is equivalent to mTi.(j)(p) = A0(p), 
with the identifications 



c(c-b) 



m = fihuj 
1 



This leads to the energy eigenvalues for ( ^.5| ) 



o = l b = 13 

A + mca = - — 

Tiuj 



(2.11) 



Eu — Tiuj 



(2.12) 



which agrees with the past results. It is a simple exercise, too, to obtain the wavefunctions. 

Note that if we wish to use shape invariance to solve ( |2.5| )when the term linear in 
is explicitly present, it is still possible to do so, as we now demonstrate. 

Suppose we define modified raising and lowering operators 



A = F(p)-^ + W(p) + n(p) 
dp 

A = -F(p)4- + w(p)-n(p) . 

dp 



(2.13) 



Then 

H = AA 

H = AA 



F{p) 



F{p) 



dp 

d 
dp 



F{p) 



+ F(p) 



fdW dn\ 



\ dp dp J 
fdW dn\ 



d 



+ — - 2F(p)n(p)— + w z ( P ) - n z ( P ) 



dp 
d 



To make contact with the problem at hand, we set 
F(p) = (a + bp 2 ) W(p) = cip 



Q(p) = C2P 



(2.14) 



(2.15) 
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in which case 



d \ d 
H = - ( (a + bp 2 )—) -2(a + bp 2 )c2p— + (c 2 -c 1 b-cl-C2b)p 2 -c 1 a-c 2 a ,(2.16) 
dp J dp 



which now includes the ^ term. Computing TC = AA, we find 

TC(a, b, ci, C2) = H(a, 6, c\ + 6, c 2 ) + 2ac + ba . 



(2.17) 



This result is identical in form to the shape invariance result (|2.9|) ; the addition of fi = c 2 p 
thus has no bearing on the eigenvalue spectrum. Even in considering the wavefunctions, 
there is no physical modification associated with including this term proportional to 7, 
simply a change of variables. 



3. Multi-dimensional Oscillator 

We now generalize our results to another shape invariant eigenvalue equation. This 
will turn out to be the relevant equation for the "radial" (in momentum space) piece of 
the generalized Schrodinger equation that describes the .D-dimensional harmonic oscillator 
with a minimal length uncertainty relation, the second of the two models discussed in H . 

To obtain our new shape invariant model, we define 



d 



9 



A = (a + bp 1 ) — + cp+- 



dp 



P 



A 



d 
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(3.1) 



a + bp )- — h cp H — 



dp p 

Then we can compare 7i = AA and 7i = AA. One finds 



H(a, 6, c,g) = - 
H(a, 6, eg) = - 



(a + bp 2 )-^- 
dp 



+ g ( 9 t a) + c(c - b)p 2 + c(g -a)+ g(c + b) 



+ 



p« 
9(9 ~ a) 



(3.2) 



P 



+ c(c + b)p 2 + c(g + a)+ g(c - b) 



The shape invariance is obvious, with 



Ti(a, 6, c, g) = TC(a, 6, c + 6, g — a) + 4a(6 + c) — 4bg . 



(3.3) 



The eigenvalues for the equation H4>k(p) = h~k4>k(j>) are readily obtained algebraically. 
They are given by 

A fc = Aabk 2 + 4(ac - bg)k . (3.4) 
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As it stands, this equation is not quite the momentum-space radial equation for the 
modified D- dimensional harmonic oscillator. That equation, given in 0, is 



-\xhuj 

p2 



[l + ((3 + (3')p 2 ] 



— V 

dp J 



+ 



D - 1 

P 



+ ( 7 D - 2(3L 2 ) + { 7 (/3Z> + + 7 ) - p 2 L 2 }p 2 



1 2E 

Rip) + ^-p 2 R{p) = t- 



d 
dp 

R(p) 



(3.5) 

which is obtained from quantizing the classical D-dimensional harmonic oscillator Hamil- 
tonian H = ^fiu 2 x ■ x + j^p-p according to the modified commutation relation 

[xi,pj] = ih(6ij + (3p 2 5ij + (3'piPj) , (3.6) 

which is done using the representations 



Xi = %h 



d 



d 



(l + (3p 2 )—+ (3'p lPj — + 7 p. 



(3.7) 



dpi dpj 
Pi = Pi 

Note that, as in the previous example, there is no physical content to 7 and, it may simply 
be set to any value. 

The problem, however, in comparing (|3.2|)with (|3.5|), is that TC in the model we have 



written down has no term linear in -j- , yet such a term cannot be removed from the multi- 
dimensional oscillator equation simply by re-defining 7, as was the case in one dimension. 
Fortunately, the same method used to include the term linear in ^ in the previous section 
can be used productively here. 

We once again define modified A and A operators, namely 



A 



a + bp 2 )-^- + (ci + c 2 )p + 



9i +92 



dp 



d 



A = -(a + bp 2 )— + (ci - c 2 )p + 
dp 



P 

9i ~ 92 
P 



(3.f 



Defining, as usual, TC = AA and TC — AA, we obtain 

2 

TC 



9i 



(a + bp 2 )— - (a + 6p 2 )(ci - ^ ) - (a + bp 2 )(c 2 
dp J p A 



-2(a + bp 2 )(c 2 p+—)4- + (ciP+—f - (c 2 p + 

p dp p P 



92, 2 



TC = - (a + bp 



d_ 

dp 



(3.9) 



+ (a + bp^)( Cl 



91 \ ,_ , l._2a/„ 92 



pz 



-2(a + bp 2 ) (c 2 p + 9 -)^ + (c lP + 9 -±) 2 - (c 2 p + 9 -l) 2 
p dp p p 
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If C2 = #2 = 0, then we simply recover the operators of ( |3.2| ); however, the terms with 
C2 and #2 in (|3.9|) are the same in both 7i and 7i, and thus there is a shape invariance 
relationship between H, and Ti in this new scenario which is identical to the relationship 
that followed from ( [3.9| ), even when C2 and g^ are non-zero. Note, though, that with c-i and 
#2 non-zero, 7i and 7i have exactly the same mathematical form as ( |3.5|) , the equation for 
the minimal uncertainty relationship .D-dimensional harmonic oscillator. Consequently, 
we can use the shape invariance of the new 7i and 7i to obtain the solutions to ( |3.5| ). In 
particular, we have 

7Y(a, 6, ci, c 2 , gi, #2) = H(a, 6, ci + 6, c 2 , #1 - a, 5(2) + 4a(6 + ci) - 46^i 

(3.10) 

A fc = 4a6/c 2 + 4(ac - 6fir)/c . 

These A& give the spectrum for the modified D-dimensional oscillator, with the appropriate 
identification of parameters. Likewise, the ground state is found by setting A<fi = 0, and 
the excited states by the application of A. 



4. Additional Shape Invariant Models in this class 

In this section, I will describe the extension of the shape invariance techniques devel- 
oped here to models that fit into the present analysis, but that go beyond that considera- 
tions of || and the simple minimal length uncertainty principle. 

Suppose we have an eigenvalue equation Hep = E<fi where H is given by 



H = - 



(a + frsinh 2 y) 



d_ 

dy 



+ 



7 



cosh y 



(4.1) 



and the parameters a, 6, and 7 are all positive. Here, I have deliberately labeled the variable 
y to emphasize that the technique described herein is not restricted to a particular physical 
framework. Let us define 



A = (a + b sinh 2 y)-f- + g tanh y 
dy 

A = —(a + b sinh 2 y)— — h g tanh y 

dy 



(4.2) 



Then 



H(a,b,g) = AA = - 



[a + b sinh 2 y) 



A. 

dy 



-1 2 



g(g + a- b) 
cosh 2 y 



+ g(g-b) 



(4.3) 



where g, like a and b, is positive. The partner operator TL = AA is related to TL by 



H{a,b,g) = H(-a,-b,g) = H{a,b,g + b-a) + (2g-a)(a-b) . (4.4) 

Thus this model exhibits shape invariance under the transformation g — > g + b — a. Using 
the standard techniques, we see that the eigenvalues Ek for the operator TL are simply 

E k = k(2g-b)(b-a)-k 2 (a-b) 2 k = 0,1,2,... (4.5) 

Note that for large enough k, these energies become negative, which we know is not phys- 
ically allowed. This is not special to the example at hand. In fact, if we consider ordinary 
quantum mechanics with potential 1/ cosh 2 (x), there is only a finite number of bound 
states; these are the states extracted by the standard shape invariance argument, and 
when the energy goes negative, it is the signal that we have left the realm of physical 
states. There is always at least one normalizable state, which is the state with vanishing 
energy in the supersymmetric/shape invariance framework. This entire structure persists 
in the problem at hand. 

Note, too, that it is simple to make contact with the original operator H in ( |4.1| ) 
above. To be precise, 

H{a,b,j) = H{a,b,g)-g(g-b) (4.6) 

where g = ~[(b - a) + a/(6 - a) 2 + 4^ ]. 

Finally, we note that if generalize this procedure to equations given by TL = BB where 

B=(l + f(y))j-+g(y) B = - (l + f(y)) j- + g{y) , (4.7) 

then we can always identify models in which 

m = (4. 8) 

for which BB has the same mathematical form as BB, and the spectra are positive semi- 
definite, so that the shape invariance approach is applicable. 



5. Conclusions 

In this paper, I have demonstrated that shape invariance can be applied effectively 
to exactly solvable models in extensions of quantum mechanics. These include models 
that possess a modified uncertainty relationship, namely the minimal length uncertainty 
relation, which also emerges in string theory. As usual, shape invariance provides a simpler 
and clearer understanding of the exact solvability of the models in question. We have 
found, too, that shape invariance can be effectively applied even for spectral problems in 
which the differential operator being studied must be written as the product of lower order 
differential operators that are not adjoints of each other. 

We note, too, that given the connection between centrally extended supersymmetry 
and shape invariance demonstrated in ||, the present work implies a supersymmetric 
construction that can incorporate these minimal-length uncertainty principle extensions of 
quantum mechanics, and the use of BPS techniques to extract the states in these theories. 
The exploration of those issues, however, lies beyond the scope of this paper. 

6. Acknowledgments 

This research was supported in part by the National Science Foundation under Grants 
No. PHY04-57048 and PHY05-51164, by the Japan Society for the Promotion of Science. 
Portions of this work were conducted at RIKEN (The Institute for Chemical and Physical 
Research) and at the Kavli Institute for Theoretical Physics, both of which I thank for 
their hospitality. 



9 



References 



[1] L. E. Gendenshtein, JETP Lett. 38 (1983) 356. 

[2] F. Cooper, A. Khare, and U. Sukhatme, Phys. Rep. 251 (1995) 267. 

[3] L.N. Chang, D. Minic, N. Okamura, and T. Takeuchi, Phys. Rev. D65 (2002) 125027. 

[4] D. J. Gross and P. F. Mende, Nucl. Phys. B303 (1988) 407; 

D. J. Gross and P. F. Mende, Phys. Lett. B197 (1987) 129; 

A. Kempf, G. Mangano, and R. B. Mann, Phys. Rev. D52 (1995) 1108 (1995). 
[5] M. Faux and D. Spector, J. Phys. A37 (2004) 10397. 



10 



